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1 Introduction 



The prototype gauge/gravity duality relates finite temperature conformally invariant = 4 
supersymmetric Yang-Mills theory to 5-dimensional AdS space (times S5) with a black hole 
[U [2]. To extend this duality to usual SU(A''c) Yang-Mills theory at finite T, on which there 
is ample lattice Monte Carlo |3l HI El |6] and even (with quarks included) experimental data, 
one of the first tasks is to see how conformal invariance can be broken, how new scales can 
be introduced. A large number of models have been suggested for this |71 El [9l [ini IH ESI [131 
[m [151 [IHl [IZ|- This article is one more, very closely related to [T2] . 

The basic building blocks of the gravity dual discussed here are {z is the 5th coordinate, the 
boundary is at z = 0) the conformal factor of the 5d metric b{z) (=C/z, C = AdS radius, in 
the conformal case), the black hole factor f{z), f{zh) = (= 1 — z^/z^ in the conformal case), 
the scalar field (j)[z) = log X[z) (= in the conformal case) and the scalar field potential V{(t)) 
in the gravity action (= VI j C? in the conformal case). Assume now the boundary theory 
one is searching gravity dual for has a running coupling g'^{fJ-) with a known beta function: 
IJ,dg'^{fi)/dfj, = P^g^). The key idea now is to firstly identify ^ (in some units) with 6(z) so 
that large energy scales, ultraviolet, UV, correspond to z — > 0, h[z) ~ l/z ^ 00 and small 
scales, infrared, IR, to large z. In the IR h may approach a constant or vanish. Secondly, one 
identifies A = e*^ with . One thus has, in addition to three Einstein equations, one more 
equation bdX/db = /3(A) relating the metric and the scalar. In the analysis of ^0j-|13' the 
scalar potential V{4') in the gravity action is so constructed that a desired beta function is 
obtained in the UV and confinement in the IR. 

In this article, we take an extreme view and use the two Einstein equations not containing 
V and the equation bdX/db = /3(A) to completely solve the problem and simply regard 
the Einstein equation containing V{4') as an equation determining V^cp). We shall mainly 
concentrate on studying a simple powerlike beta function ansatz, /?(A) = — /3oA'^. Here q is 
a parameter, q > 1 corresponds to a theory confining in the IR, q = 2 would correspond to 
one-loop Yang-Mills running and q = 10/3 has in [18j been shown to reproduce the main 
features of the detailed analysis in [13 . Further, we shall modify this ansatz so that its large 
A behavior is constrained to be — | A(l -|- a/ log A -|- ..). First order transitions are obtained 
for q > 1 and a > 0, a = describes a continuous transition. 

The theoretical basis of this model has been laid in ^10j-[13j, starting from a potential V{(p)- 
The chief virtue of the approach in this article, starting from /3(A), is that semianalytic formu- 
las, requiring only numerical evaluation of integrals, are obtained. Details of SU(A''c) gauge 
theory thermodynamics can be described with q, a and two more parameters: a dimension- 
less C? jG'o fixed by p{T)/T'^ at large T and an energy scale A fixed by Tc, the transition 
temperature between high and low T phases. Note that, in particular, the numerical value 
of the parameter (5q in the beta function ansatz, never enters for physical quantities. 
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2 The model 



One starts from the gravity + scalar action (in the Einstein frame and in standard notation) 



5 = ^^11 d'x^[R-l{d^cl,f + V{<P)] -2 J d'xV^Ky 



By writing 

it can also be written in the string frame. One now assumes a metric ansatz 



ds^ = h\z) 



-f{z)dt^ + dx^ + 



(1) 
(2) 

(3) 



The four functions b{z),f{z) in the metric, the scalar field (j){z) and the potential V{(j){z)) 
are then determined as the solutions of the three field equations following from ([T]) : 



6^ b bfb'^,. 



4+3- = 0, 
/ b 

{b = b'{z), etc.) and from a fourth equation, 



dX 



(4) 
(5) 
(6) 



/3(A) = 6—, A(z) = e 



Hz) 



(7) 



where /3(A) is the beta function of the field theory one is seeking the gravity dual for. To 
begin with, we use the very simple beta function 



/?(A) = -/JoA" 



(8) 



and solve the equations in the order Q, ([s]) and ([6|. After this, Eq. ^ gives the potential 
V{(p). In Sections |5] and |6] slightly modified beta functions are studied. 
In |12| [13] the starting point was the potential 



V{<t>) 



g{l + FoA + yiA^/^[log(l + y3A^)]V2} 

where the parameter values used were 

Vq = 0.04128, Vi = 14.3, V3 = 170.4, A(z = 1) = 0.0242254. 



(9) 
(10) 

(11) 
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Vo matters only in far UV and can be set to zero for thermodynamics, which is dominated by 
the term [18]. The equations of motion (j4|-(j6| were solved numerically in |13j and (j8| can 
then be used to determine the beta function |18) . Here we proceed in the opposite direction, 
starting from the beta function. The prototype values of q would be g = 2 corresponding to 
the Vo term in (10) and q = 10/3 corresponding to the thermodynamically dominant A'''''^ 
term (see Eq. (29) below). In general, q> 1 and the limit q = l plays a special role, e.g., the 
latent heat vanishes in this limit. 



2.1 The fourth equation. 

First, from Eq.([7| and for q> \ 



log 



6o {q - l)/3oA«-i 



(12) 



where 6o is a constant, the analogue of Aqcd; the scale at which A diverges. It will prove to 
be convenient to use Q as the extra dimensional variable and relate 6, A, / and z to it. The 



numerical value of /?o and the normalisation of A do not matter, only the combination in (12) 
enters. We shall also abbreviate 



16 



a = 



% - 1) 



2 



3(<?-l) 



(13) 



An important characteristic of (12) is that in the IR, A oo, Q ^ 0, 6 approaches a non- 



zero constant 6o > 0. We shall in Section [5] study the case when h vanishes proportionally to 
a power of Q and in Section [6] the case of h vanishing ~ powers of Q and logQ. Then the 



simple relation (12) between h and Q is replaced by the more complicated one in (73) 



2.2 The second equation. 

To integrate the second equation one introduces 

W = -h/b^, 

which simplifies ([5| to the form 



(14) 



hW 



4 \2 



AVA^. 



Writing W = dW/dX ■ A cancels one power of A and replacing the remaining A by the beta 



function using the basic equation ([7| just produces to the right hand side the definition (14) 
of W. The equation then integrates to give 



l)2log(6/6o) 



(15) 



W{X) = W{0) exp 

The normalisation iy(0) 
ically AdS, V{0) = VI j C? (see Eq. (27) below). Now that W is known, the second integration 
leading to 6 = 6(z) is performed by writing ( 14 ) in the form 



1 /£ follows from the requirement that the boundary be asymptot- 



dz 



db 



dQ 
-bW 



C 

~dQ exp 

00 



-Q 



a 

4Q 



(16) 



3 




Figure 1: Plot of the integral e'^I{Q,a) for two values of a. 



or, if one wants A = X{z), 
Here the important energy scale 



dX _db(3 
dz dz b 



A 



-P{X)b{X)W{X). 

b_o 
" C 



has appeared and one has 
where 

I{Q,a) = I dyexp 



Az = I{QAa 



4 



r 

JQ 

oo ^ 

5^-(-arr(l-n,Q) 





exp 



-Q 



Q 



a 2a + 6a 6a^ + 24a 



/(0,a) = 2^Ki{2./^), 



(17) 

(18) 
(19) 

(20) 

(21) 

(22) 
(23) 



A few terms of the asymptotic expansion at large Q are exhibited here. The behavior of the 
integral I{Q, a) is as shown in Fig. [T] The Q = limit together with (23) implies that 



(24) 



a "wall" has been dynamically generated via the assumption ([7]). 
2.3 The third equation. 

The third equation (together with the first one) always has the trivial solution / = 0, / = 1. 
This will be the low temperature phase with zero free energy. It is confining for q > 1: 



4 



0.8 



Figure 2: Field configurations for q ~ 10/3 as functions of the extra dimensional coordinate. First 
Q — Q{z) is computed from (19 1, then b{z) is obtained from (12 1, z extends up to ( p4| . Then 
A = A(z) = l/[{q~'l)(3QQ{z)]^/'^^ is obtained from Q{z) by inverting ^ and setting {q-l)Po = 0.3. 
These give the string frame bs{z). Finally f{z) follows from (26 1 setting Qf^ = 1. The choices of 
parameters have no effect on the thermodynamics. 



the confinement criterion is that in the string frame the metric factor bs = b\^^^ have a 
minimum at some z (see Fig. |2]). Using ([s]) this converts to the condition that the equation 
f3{X) + I A = have a solution. This is so for q > 1 and one finds ^8], e.g., the string tension 



^~ 27ra' \2(3o) 



(25) 



To have a nontrivial solution, one changes dz to dQ using (16) and directly integrates the 
third equation ([6| to give 



C1+C2 



dz 



Ci + C2^I{4Q,a) 



1 



dz 

w. 



I 



dz 



1 



/(4Q,a) 
I{4Qh,a) 



(26) 



where C2 is fixed by introducing a scale Qh = Q{zh) at which / vanishes, f{Qh) = 0, and 
Ci = 1 by the fact that b{z) = C/z, z ^ demands f{z = 0) = f{Q = 00) = 1. 
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1/e = iq-DPo^"- 



0.0 



2.0 



Figure 3: V{Q) from (28 1 for q = 10/3. The dashed hne is the full potential from [13 , indistinguish- 
able from the Q > Qh limit ( 29 1 . 



2.4 The first equation, tlie potential V{(f)). 

Now that we have all the functions solved, we can insert them to the first equation Q and 
obtain the potential 



l2fW 
12 / a \ 



12 
£2 



1 + 



+ 



3A 



1 - 
2a2 



f 
b 

/(4Q,«) \ 

I{AQn,a))^''^\AQ 



8Q 16Q2 



+ 



12 
£2 



1 + 



9{q-l) 



+ 



exp(-4Q) 
/(4Q/^,a) 



(27) 

(28) 
(29) 



Eq. (|28[ is plotted in Fig. g for a = 16/(9(g - 1)^) = 16/49 by choosing some values 
for Qh- V is always positive in the form (28): there is a range beyond the horizon where 
/(I — (/3/(3A))2) < but the last term cancels this negative part. One sees how the potential 
agrees with the potential in |T2l [13] for Q > Qh, z < z^. Analytically, the Q ^ Qh, 
A ^ 1, z <^ Zh limit (29) coincides with the same limit of the potential in |12| IT3] . Beyond 
the horizon V starts growing rapidly, in fact, for Q <^ Qh, V{Q) ~ Q~'^ exjp[a / {2Q)] ~ 
A^-^-^expig^^/JoA^-i]. 

The solution b{z), (j){z), f{z) for the thermal case is given concretely by the equations 
above. To get the vacuum solutions bQ{z), 4>q{z), / = 1 one should now solve Eqs. (Q and 
^ with / = 1 and using the potential (p8|). We do not need these solutions, only that they 
exist. 



3 Thermodynamics 



Now that f{z) is explicitly known in (26), the standard formula 47rT = —f'{zh) with dz/dQ 



from (16) gives the horizon temperature (when Qh is a dummy variable we leave out the 
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index h). Thus 



AttT 



dQ 



or 



7rr(g) 

A 



exp(-3(5) 
I(4Q,a) 



+ 



3a 



16Q2 32Q3 i28Q^ 



12a + 
^ 512Q5 ^ 



(30) 

(31) 
(32) 



Conversely, Eq. (31) gives Qh = Q?t('/rT/A). The large-T asymptotic expansion of this 
relation is, inverting (32) and denoting L = log(7rT/A), 



Qh = L 



a a 
4L 16L2 



a + 2a2 12a + 25a2 

+ — + ... 



32L3 



512L4 



(33) 



A particular role is also played by dT/dQ. In fact, below in (43) we shall prove that 
dT/dQ = 3cl{T)T. 

The relation T = T{Q) is plotted in Fig. |4]in the small Q region. At large Q the asymptotic 
expansion (32) is very accurate and T ~ exp{Q). At Q = 



7rT(0) 



vr dT{0) 



A 2^i^i(2^)' A dQ 



-3^T(0). 



(34) 



Thus dT/dQ has to change sign at some Q = Qmin. At this point there is a minimum of 
T{Q) and = 0. In the more general analysis of [12] the branch with dT/dQ < is called 
the small black hole branch, the one with dT/dQ > is the big black hole one. In Section 
|6] we shall discuss the case when T can grow to infinity in the small BH branch, now it is 
bounded by ( 34 ) . 



From the form ([3]) of the metric it follows that the black hole entropy density is 



(35) 



Scaling with T, from (31), (32) and (33), 
(7r£)3 



s 

f3 



AG, 
AG, 

AG, 



[exp(4Q;,)/(4Q;„a)]^ 



exp 



1 



—3a 
4Ql 



1 + 



3a 



3a 



IQQi 32Qi 



+ 



3a^ + 9a 15a^ + 36a 



256Q^ 



h 



512QI 



+ 



3a 9a2 + 6a 9a^ + 2Aa^ + 12a 
iZ ^ 32L2 128L3 ^ 



L = loe 



vrT 



(36) 
(37) 



The full bulk thermodynamics is known when the (minus) free energy density or pressure 
p{T) is known. The rest then follows from relations like 



s{T)=p'{T), e{T) = Ts-p, 



3p 



2^4 



T — c 

dTT^' ' 



dp 
~de 



Ts'{T) ' 



(38) 
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Figure 4: T{Q) from (31 1 (left panel) and p(T)/T'* from (39 1 (right panel) with numbers correspond- 
ing to q = 10/3, a — 16/49. The pressure vanishes at p{Qc = 0.363) = 0, the corresponding value of 
T is ttTc/A = 1.67. At this point = 0.707/3. The dashed line is the leading large Q approximation 
exp[Q + a/(4Q)]. 



Now that s{T) is known, p{T) can be obtained by integrating s{T) = p'{T). The key 
observation in [12|[T3] now is that the correct constant of integration is obtained by integrating 
over the entire curve T = T{Q) starting the integration at the value T(0) corresponding to 
Q = in Fig. |4} Thus one has 



p{T) 



T 



A3 rQ{T) jry 



£A. 



(39) 



Since dT/dQ starts negative, p{T) starts decreasing (see Fig. |4]), but soon starts growing 
and crosses the value zero at some T which defines the transition temperature between the 
phase we are discussing and the phase corresponding to / = 1 in the metric ansatz ([s]). Since 
this low T phase has p = 0, the transition temperature is determined by 



piQc) = piQiTc)) = 0, 



(40) 



where the value of Tc = number x A is obtained from (31 ). 



From (39) one further computes 



e{T) 



4G, 



[rexp(3Q)-p] 



e{T)-3p{T) 



4G5 
"o 

4G5 



QiT) 







dQ T{Q) exp(3Q) + T(0) 

dT 
dQ 



3 / dQ It- — ] exp(3Q) + r(0) 



(41) 
(42) 



and for the sound velocity 



1 dT 
3f dQ 



4—— exp 
A 



-Q 



a 

4Q 



(43) 
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(e -3p)lT^ 



TIT, 



Log[r/r,] 



Figure 5: Left panel: e/T*, Z'pjT'^ and the scaled interaction measure up to SOTc for q = 10/3. Here 
and in other bulk thermo figures quantity /A^^ is plotted. The dashed line shows the limit at T ^ oo. 
Right panel: p{T) at very large T scaled by the T ^ oo limit. The two data points at T = SOOTc 
and r = 3 • lO^Tc are from [2D , the third is the highest from [3 . The dashed lines are the asymptotic 



expansions (44l including only the — 3a/41og(T) term. 



The large Q or large L = log(7rT/A) expansions are 

p (7r£)3 /I 3a Oa^ Oa^ + 12a 27a^ + OGa^ 

_l 1 1_ 



T4 4G5 V4 16Q 128Q2 512Q3 8192Q4 

(7r/:)3 ^1 3a 9a2 9a3 + 24a2 + 12a \ 



4G5 V4 16L 128L2 512L3 
e _ (-nCf fl 3a 9a2 + 8a 9a3 + 240^ + 12a 
r4 ~ 4G5 U ~ 16Q ^ 128Q2 512Q3 ^ 

(7r£)3 /I 3a 9a2 + 8a 9a3 + 48a2 + 12a 

U ~ 16L + ^28L2 512L3 ^ " ' ' ^ ^ 

e-3p _ (7r£)3 / a 3a2 90^ + 80^ + 12a 

~ 4G7 V16Q2 ~ 64g3 ^ 512Q3 ^' 
(7r£)3 / g 3a2 Qa^ + 24a2 + 12a \ 
4G5 64L3+ 512L4 +7' 

" T7(r) "3V 4Q2+^ + ■■■) ■ ^^^> 
Note that to plot the asymptotic expansions vs T/Tc one must write 



and determine vrTc/A from (40). 
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012345012345 



Figure 6: e/T"', 'ip/T'^ and the interaction measure up to ST^ for q = 10/3 (left) and q = 2 (right). 
For g = 10/3 the interaction measure decreases ~ 1/T^-^ above Tj, (dashed line), for g = 2 it decreases 
~ 1/T^ for Tc < T<2Tc (dashed line), in agreement with data, but more slowly for T>2Tc (see also 
Fig.ll. 



0.6 r 




2.0 2.5 3.0 3.5 1.5 2.0 2.5 3.0 



Figure 7: Left panel: The latent heat plotted as a function of q. L vanishes at g = 1. The lattice 
value at iVc —> oo is 0.34(6) [4J and is obtained for q = 2.2. Right panel: Lattice data for Nc ^ I 
(dashed lines) [6] compared with model prediction for q = 2.2. 

4 Comparison with lattice data for hot SU(A^c) gauge theory 

Our model is so far based on a very simple powerlike in A (exponential in (p) ansatz for the 
beta function, from which the scalar potential is derived, and one cannot expect detailed 
agreement with data. Nevertheless, we find that on a qualitative level the main features of 
the data are reproduced to surprising extent. 

To analyse the bulk quantities one has to fix the value of the parameter /G^. We fix it 
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Figure 8: Left panel: Data and model with q — 2.2 for the interaction measure scaled by . Right 
panel: "i-pjT^ ioY q — 2 plotted vs. TjT^. Continuous curve (exact), dotted curve (Eq. (44 1 up to 
1/log), dashed curve (Eq. (44 1 up to 1/log^). 



to the pressure of ideal gluon gas, 



PSB = 2(iV2-l)^r^ (49) 



Comparing with (44) one has 

This is less by a factor 8/45 relative to the value for = 4 SYM, 2/15 comes from the 
reduction of the number of degrees of freedom and 4/3 from going from strongly interacting 
to ideal systems. Bulk thermodynamic quantities are very accurately proportional to A'^^ 
in [5l |6] so that we shall simply divide by this factor. 

Using this normalisation Figs, [s] and [g] show e/T^, p/T^, (e — 3p)/r^ as functions of T/Tc 
evaluated numerically from the above equations for q = 10/3 and q = 2. Fig. [5] also shows 



p{T) at very large T scaled by the T ^ oo limit (49). 



At Tc, p = p(Tc) = but e(Tc) has the nonzero value given by (41). Since p = e = for 



the low T phase, p is continuous, by construction, but e jumps. This gives the latent heat 
^ - ^ = ^ = '-^N^ [exp(4Q.)/(4Q„ a)f . (51) 

This is plotted in Fig. [7] as a function of q. The latent heat vanishes for ^ = 1, the transition 
goes over to a second order one. The lattice value [1] L/T^ = 0M{6)N^ would be obtained 
for q ~ 2.2. Choosing this value of q, Fig. [7] also compares the T dependence of the bulk 
quantities with lattice data extrapolated to A"c — > oo [6j. 

One well known characteristic of (e — 3p)/T'^ is that it decreases ~ 1/T^ for T>Tc. This is 
checked in Fig. plfor the data and for the model with q = 2.2. As seen in Fig. pi for q = 10/3 
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the model result decreases faster, ~ 1/T^-®. Continuing to still larger T, according to (46) 
the leading large T behaviour of the interaction measure is 

e - 2,p 4tt'^N^ 1 



r4 135(g - 1)2 log2(^r/A) ' 
a slight generalisation of Eq. (1.3) in |12] . 



(52) 



The asymptotic expansion (44) looks much like a perturbative expansion in and it is 
of interest to compare the leading corrections quantitatively - to higher orders the finite T 
perturbative expansion proceeds in powers of g. Consider first the interaction measure. In 
perturbation theory, to leading order, by taking Td/dT of Eq. (6.1) of |19j . 

e-3p_ {N^ - l)N,^ dg\T) _ TT^jN^ - 1) 1 ^^^^ 



144 dT 66 log\7TT/A^y 



This is very much similar to (52), in fact, coincides for q = y^88/45 + 1 ~ 2.4. Consider then 



the pressure, for which the two expansions are 

%^ =1 4 

PSB 16^2 ^\ I67r2 J ^ 3{q - 1)2 log(7rr/A) ^ ' 

where also the g^ plasmon term is included in the perturbative expansion. Comparing the 
leading terms with 1-loop running for (?^(T), the two expansions again coincide for q = 
Y^88/45 + 1 ~ 2.4. However, the plasmon term is quantitatively very important, smaller 
than the g^ term only for T > lO^Tc, which explains why the q = 10/3 term agrees with 
the perturbative result, see Fig. |5j This figure also shows to what extent the leading order 



term of the expansion (44) agrees with the numerical computation. How still further terms 
improve the agreement is seen from Fig. |8j 

From the above one sees that q = 10/3 gives a good fit to SU(A'^c) thermodynamics for 
T>10Tc; for smaller T the curves are too high and the latent heat too big. This indicates 
that one has to modify the ansatz in the IR large A region, the excellent fit in [13] contained 
several parameters. In the next two sections we present systematic generalisations of ([s]), 
which will be seen to lead first to vanishing latent heat, a continuous transition and then 
again to a first order transition. 

5 Continuous or second order transition 

As discussed above, the transition approaches a continuous one when q ^ 1. To study this 
case it is more realistic not to take g — > 1 in ([8| but to take a beta function which only 
approaches the marginally confining q = 1 form — | A at large A: 

Even for this case the solution of the gravity equations (|4])-([7]) is surprisingly simple. Denoting 
again (but now this is not log(6/6o)) 
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and recalling the abbreviations 



A 



we have 



W 



60 
1 



3(g-l) 



Q_ 
Qo 
A 



1 + 



Q 



(57) 

(58) 
(59) 



where Qq is the scale at which b has the value ftoe*^"- In Section [g] we shall choose Qq = A. 
Note that now b vanishes in the IR proportionally to a power of Q or of 1/A. 

Using db/dz = —b'^W from (16) and A-kT = —f{zh) from (26) one finds the temperature 



A 1 



JiQ q3A 



-3A-1 



1 

fiQ,Ay 



(60) 



the pressure (Qo cancels from here) 



p_ 



(47r£)3 J_ 

4G5 r4 



dQ 



dT 



and sound velocity 



1 Q 



dQ 
dT 



,3Q q3A 



(61) 



(62) 



' 3 A + Q TdQ 

A numerical example of the scaled bulk quantities for q = 10/3 is shown in Fig. |9j 

Now that the transition is of second order the minimum of T(Q) is at Q = 0; beyond that 
dT /dQ > 0. T{Q = 0) is calculable and defines the critical temperature 



ttTc 
A 



A 
Qo 



(63) 



To find the behavior of T{Q) near Q = or near Tc one must do partial integrations with 
dy~^^ to increase the power of y in ( 60 ) . For < 3^ < 1, q > 3 one partial integration is 
enough to make the integral converge for Q — > 0. Expanding the remainder in Q one obtains 

f{Q,A)/f{0,A)-l 
= A^-' 



dyy 



-3A g-4j/ 



T/Tc 
5A 



1 = t 
l + 4y 



{A + y)A 



Q 



3A 



6A 



1-3A 



Q+om, 



(64) 



where the integral is positive for < 3^ < 1 . For this range of 3 A the first term is dominant 
in the critical region t — > so that 



tr^Q 



3A 



Q^t 



1/3A 



< 3A < 1. 



(65) 
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For 1 < 3 A < 2, 2 < q < 3 two partial integrations are needed and 



t = ^^—^ Q - Integral Q^^ + O(Q') , (66) 



where the Integral, similar to the one in (64), is positive. Thus now and for all larger values 
of A and smaller values of q: 

tr^Q, 1<3^<2, 2<g<3. (67) 

When converted to the critical behavior of the free energy (= —p) one has 

3A<1 f^t^^+\ 3A>1 (68) 

so that the critical exponent in Cy ~ f"{i) ~ is a = (= 1 — 3^) for 3^4 < 1, {3A > 
1). Since f"{t) is not divergent, the transition is a continuous one. For the SU(2) finite 
T transition, of 3d Ising universality class, one has a = 0.12 [21], a mild divergence. A 
continuous transition appears also if one takes SU(3) gauge theory with effectively 3 flavours 
of infinite mass fermions and reduces the mass of one fermion, i.e., studies one-flavor QCD. 
For sufficiently small mass the transition becomes a continuous one f22^ , the end point again 
being of Ising universality class. 



With the modified beta function (55) with q = 10/3 we have now reduced the latent heat 
) zero. The work in \TI\ [TB] 
:der transition but still mair 
behavior at large A should be 



to zero. The work in \TI\ [TB] indicates what the modification should be if one wants a first 

3 
2 



order transition but still maintains the leading asymptotic behavior — |A: the asymptotic 



+ (69) 

witt[^ a > 0. We shall in the next section derive analytic formulas even for this case. In 
principle, one could then redo the critical analysis above so that the parameter a would 
effectively be an external field which could drive the transition to a first order one and more 
critical indices would appear. 



6 First order transition again 

Consider now the following ansatz: 



/3(A) 



1 + 



A 1 



-1 



q-l Q+A 



l + a( 



l + a( 



1) 



^ log(l + I /3oA'?-i) 
'log2(l + |/3oA''-i) + l 

log(l + A/Q) 
'log\l + A/Q) + l 



(70) 



(71) 



where Q, A are defined as in the previous section. The modification factor relative to ( |55[ ) is 
so constructed that in the IR, A — > oo, it is 1 + a/ log A + in the UV it is 1 + 0{\'^^^) and 



also so that W is integrable in closed form from (15): 

A 



w 



1 



1 + 



a 



(72) 



'^In the notation of T?, this parameter was denoted by | a/{a — 1), a > 1 
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Figure 9: Scaled bulk thermodynamic quantities for the beta function (56 1, leading to a continuous 
transition, with q — 10/3. If the curves were normalised as in (50 1, they should be multiplied by 
2567rVl5 « 168. 



The expression for h from ([7|, 

^ b d\ dQ -A , A 1 



^ log2(l+A/y)+l 

is also surprisingly simple, though not integrable in closed form. The lower limit of the Q 
integral has here been fixed to be A, i.e., h{Q = A) = bo. In the IR, Q ^ 0, this implies 

This is the key modification relative to ( |58| ) : b vanishes in the IR but more slowly by a power 
of logQ than an ansatz with /3 — > — | A. 

The modified IR behavior has several consequences. Consider first the z dependence of the 
field configurations. For this one has to integrate z = z{Q) from dz = —db/{b'^W) so that, as 
a generalisation of (19), 

Jq{z) dQ bW 
Studying the limit Q{z) —>■ one sees that this integral diverges if | a < 1: 

Az~log(^/Q)i- t", (76) 
otherwise it is finite. Thus, if a < | the range of z is unbounded. 
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Similarly, for the temperature we have 



and counting IR logs one finds that the power of log(l/(5) is 2a from the outside 6^, 
-8a/3 - 2a/3 from the inside l/{b'^W) so that in the IR T ~ log''"/^(l/Q). Thus T di- 
verges logarithmically when Q ^ 0; with the beta function (|8| T{Q = 0) was finite but 
there was a minimum, Fig. |4j Thus there again is a minimum of T{Q) and a branch with 
dT/dQ < 0, so that p{T) behaves as in Fig. |4]and the transition becomes of first order. Full 
thermodynamics in terms of two parameters, q > 1 and a > 0, is now given by the formulas. 



generalisations of (39) and HI 



s{T) = ±-b\QiT)), (78) 
p(r) = j ^Tsm=-l ,Q-m), (79) 

<T) = — dQT{Q)b'{Q)—^, (80) 



4^5 70 dQ 
3 /"^^^^ fr^dlogb dT\ 3 



2 1 dT dQ 1 / ttT \ 



Note that r(0) in (41) is actually T(0)b^{0) which now vanishes. The procedure again is that 



one first finds at what value of Q = Qc the pressure vanishes: p{Qc) = 0. From (77) one 
next finds the corresponding value of T = Tc = T{Qc). Finally one plots parametrically with 
Q > Qc SuS the parameter the required thermodynamic quantity as the y axis and T{Q)/Tc 
as X axis. 

We refer a detailed comparison with lattice data in the few range and with QCD pertur- 



bation theory at large T to later work. As a preliminary step. Fig. 10 compares the entropy 
density scaled by N^T^ obtained for q = 10/3 with the Nc ^ oo limit of lattice data (known 
jS] for T < 3.4Tc) and the QCD perturbative result at large T. One sees that the too 
large latent heat observed in Fig. |6] for this value of q (which leads to good large T behavior) 
is cured by the modification of the IR region incorporated in the beta function (71), if one 
chooses a to be in the range 0.1... 0.2. There is still some deviation very close to Tc, but 
overall the agreement is excellent, as already noted in |13j . 

Finally, one may note that the latent heat vanishes when a — > as a power of a: 

^ = 1.55ai-=^o, q=f. (83) 



This is accurate up to a = 0.1. 

■^The perturbative result is obtained by computing entropy density from the pressure in Eq. (6.4) of |19| . 
choosing the renormalisation scale fi = 2nT, the unknown 4-loop Linde coefficient Qc — —3800 and a 2-loop 
running coupling with Nc = 3. 
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3 



TIT, 

4 



Log(r/r,) 

3 4 5 



Figure 10: Entropy density scaled by T^N^ (tlie asymptotic T ^ oo value is for the beta 

function ^ with q = 10/3 and a = 0.25, 0.125 and for T < bT^ (left panel) or T < lO^T^ (right 
panel) . 



The dashed curve gives the Nc 
perturbative result (see text). 



oo limit of lattice data [6J and the dotdashed one a QCD 



7 Conclusions 

In this paper, we have studied a simple gauge/gravity duality model for gauge theory ther- 
modynamics. Basically, the model is a modification of the thorough work carried out in 
[TUl (TTJ [121 US] : instead of starting from a scalar field potential F(0) of the gravity side we 
start from the beta function on the gauge theory side. Its chief virtue is simplicity, semi- 
analytic formulas are obtained, no numerical solutions of Einstein's equations are needed 
and the only essential parameters are the power q > 1 of A = e*^ in the beta function 
and the parameter a > associated with the approach to the infrared limit. Notable is 
also that the normalisation of A never enters, everything is expressed using the combination 
Q = l/((g-l)/3oA^-i). 

The key element in this analysis was that the three equations determining the three func- 
tions b{z), (p^z), A(z) were taken to be the two Einstein equations not containing the potential 
V{(j)) and the equation relating the beta function to b{z), A(z). The Einstein equation con- 
taining the potential then serves simply to determine the potential. 

It is perhaps surprising that such a simple beta function as (3 = — /JqA"^ reproduces main 
features of a first order transition so well. Detailed agreement with lattice data and QCD 
perturbation theory requires beta functions constrained to approach — | A(l -|- a/ log A) at 
large A. Here a > leads to a first order transition. Even for this case analytic expressions 
could be given. They become particularly simple for the case a = 0, corresponding to a 
continuous transition. 

There are many directions into which one could develop this work. In particular, a study 
of beta functions of the type — /3oA'^(l — A/A*) containing an infrared fixed point A=k suggests 
itself. Further, the model could be applied to walking technicolor beta functions which only 
approach a fixed point A* . Perhaps then one looses the virtue of simplicity. 
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